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ABSTRAGT 

We present the first results from the SAGE-Var pr ogram, a follow on to the Spitzer legacy program 
Surveying the Agents of Galaxy Evolution (SAGE; |Meixner et al.||2006|. We obtained 4 epochs of 
photometry at 3.6 & 4.5 ^m covering the bar of the Large iVIageila!nic Cloud (LMG) and the central 
region of the Small Magellanic Gloud (SMG) in order to probe the variability of extremely red sources 
missed by variability surveys conducted at shorter wavelengths, and to provide additional epochs 
of observation for known variables. Our 6 total epochs of observations allow us to probe infrared 
variability on 15 different timescales ranging from ~20 days to ^5 years. Out of a full catalog of 
1717554 (LMG) and 457 760 (SMG) objects, we find 10 (LMG) and 6 (SMG) large amplitude AGB 
variables without optically measured variability owing to circumstellar dust obscuration. The catalog 
also contains multiple observations of known AGB variables, type I and II Gepheids, eclipsing variables, 
R GrB stars and young stellar objects which will be discussed in following papers. Here we present 
infrared Period-Luminosity (PL) relations for classical Gepheids in the Magellanic Glouds, as well as 
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improved PL relationships for AGB stars pulsating in the fundamental mode using mean magnitudes 
constructed from 6 epochs of observations. 

Subject headings: 


1. INTRODUCTION 

The study of stellar variability has had a wide-ranging 
impact on astronomy and cosmology. Recent variable 
star surveys such as the MAssive Com pact Halo Object 
search (MACHO; Alcock et al. 1997) and the Optical 
Gravitat ional Tensing Experiment (OGLE; Udalski et 
^|1997|) have generated catalogs of tens of thousands 
01 variable stars of numerous classes in the Magellanic 
Clouds. However, like most ground-based variability 
surveys, both of these were performed at visible wave¬ 
lengths, and therefore miss the reddest variable sources, 
such as dust-enshrouded, highly evolved Asymptotic Gi¬ 
ant Branch (AGB) stars, which are nearly invisible ex¬ 
cept in the infrared (IR). AGB stars are unstable, and 
exhibit variability on t imescales of hundreds of days (|Vas- 


siliadis & Wood 


1993). Due to this, they are also classi- 


hed as Long Period Variables (LPVs). There are several 
surveys that monitored LPVs at near-IR wave lengths in 
the Magellanic Clouds (J, H , and AT-bands; Ita et al. 
2002 Whitelock et al. 20031), but even these can still 
miss the dustiest soured The pulsation properties of 
the heavily enshrouded (and most evolved) AGB stars 
are therefore not well understood. This problem can be 
addressed with IR monitoring at wavelengths longer than 
3 /rm, but there are very few e xamples of mid-IR mon- 
itoring surveys in any galaxy. Le Bertre (1992 19931 
obtained light curves at 1-20 fim ot ^60 O- and G-rich 
AGB stars in our own galaxy and found that while pulsa¬ 
tion amplitudes generally decrease into the mid-IR, cir- 
cumstellar dust can c ause amplitudes to inc rease again 
at A > 3^m. In M33, McQuinn et al. (2007) obtained 5 
epochs of imaging at 3.b, 4.5, 5.8, and S/iin and found 
that the pulsation amplitude tends to increase with color, 
but that this relationship may break down at the longest 
wavelengths. In the Magellanic Clouds and other nearby 
dwarf gala xies, only 2-3 epochs of imaging is available at 
A > 3 Mm (iP olsdofer et al.|2015[|Viih et al. 2009[[Boyer et 
al.|2015a|bV ' While these surveys can detect tlux changes 
m a large traction of the dustiest stars, they cannot place 
stringent limits on the pulsation periods or amplitudes in 
the infrared. In order to explore the variability of these 
reddest sources, which dominate the mass return to the 


interstel lar medium (ISM) from evolved stars (Riebel et, 
al.||2012|) , we used the Spitzer Space Telescope to survey 
the bar regions of the LMC and SMC at 3.6 and 4.5 /im. 
SAGE-Var represents the first large scale variability sur¬ 
vey at such red wavelengths. While our original focus 
was the reddest AGB stars, we have detected over 2,700 
IR variables in the Clouds of many classes. 

This paper is organi zed a s follows: In jj[2] we d etail the 
observational strategy (13, data reduction \2.2 \ , cat alog 
construction (2.3), and source classification scheme (2.4) 
for the SAGE-Var program. In §[^ we discuss our results; 
the identification of Long-Perioa Variab le (L PV) candi¬ 
dates without prior observed variability (3.1), an investi¬ 
gation of the connection between variability amplitude 
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TABLE 1 

Observation dates 


Epoch 

Date 

Julian Date 

Epoch 1 

LMC: 

2005 Jul 20 

2453572 

Epoch 2 

2005 Oct 28 

2453672 

Epoch 3 

2010 Aug 17 

2455426 

Epoch 4 

2010 Sep 10 

2455450 

Epoch 5 

2010 Dec 25 

2455556 

Epoch 6 

2011 Apr 27 

2455679 

Epoch 1 

SMC: 

2008 Jun 15 

2454633 

Epoch 2 

2008 Sep 19 

2454729 

Epoch 3 

2010 Aug 17 

2455426 

Epoch 4 

2010 Sep 12 

2455452 

Epoch 5 

2010 Dec 24 

2455555 

Epoch 6 

2011 Jun 16 

2455729 


Note. — Dates of the 6 epochs of 
observations of the SAGE-Var project, 
including the 2 epochs of the original 
SAGE-LMC and SAGE-SMC programs 
(Epochs 1 & 2). The dates listed for 
the original two epochs of observations 
are the approximate midpoints of the 
~5 day observation periods, while later 
epochs were executed during a single 24- 
hour period. 


and Dust Production Rate (DPR) from evolved stars 
(3.2), and new measurements of the IR Period Lumi¬ 
nosity relationship for bot h LP Vs and classical Cepheids 
in the Magellanic Clouds (3.3). Our conclusions are pre¬ 
sented in §1^ 


2. THE DATA 
2.1. Observations 

The 4 epochs of SAGE-Var observations {Spitzer PID 
70020) were taken over a 10 month Mriod, between Au¬ 
gust 2010 and June 2011 (see Tabled. Each epoch is a 
3.7° X 1.5° (LMC, Figure 3 or 1.7=^ 1.7° (SMC, Fig¬ 
ure 3 mapping of the bar region of the galaxy, using 
both the [3.6] and [4.5] bands of the Spitzer warm mis¬ 
sion. All frames were taken as 12 second exposures using 
the IRAC High Dynamic Range (HDR) mode, which also 
produces short 0.6 second exposures in order to mitigate 
the effects of saturation for the brightest sources. This 
exposure mode was chosen to provide unif ormity with 
the prev ious two epochs of the SAGE-LMC ( Meixner et 
ar]|2006 ) and SAGE-SMC (Gordon et al. 20lT) surveys. 
The scheduling of our observa tions was defermined using 


Madore & Freedman (2005) power-law cadence with an 


index ot 0.99, taking into account the original two SAGE 
epochs. Including these two original epochs, our 6 total 
epochs of observations allo w us to probe variability on 
15 different timescales (see §2.3). 


2.2. Data Reduction 

The SAGE-Var data were processed using a pipeline 
developed at the University of Wis consin for the 
GLIMPSE survey (Benjamin et al.|2003), which was also 
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Fig. 1.— Three-color image of the LMC (red: 8.0 fim, green: 
5.8/im, blue: 3.6/i.m) showing the footprint of the SAGE-Var ob¬ 
servations in green, focused on the stellar bar region of the galaxy. 



Fig. 2.— Three-color image of the SMC (red: 24pm, green: 
4.5 pm, blue: 3.6 pm) showing the footprint of the SAGE-Var ob¬ 
servations in green, focused on the main stellar locus of the galaxy. 


used for processing the original SAGE-LMC and SAGE- 
SMC observations. The Wisconsin pipeline corrects for 
numerous observational artifacts such as stray light, col¬ 
umn pulldown, banding, and bad p ixels. These obser- 
vational artifacts are discussed by Hora et al. (20041 
an d the Wisconsin pipelin e is described in sectioii 3.1.2 
of Meixner et al. (|2006|). The pipeline is based on 
DAUEHUT (jstetson 19871 and is modified to account 
for the variable diffuse background in IRAC images. 
Photometry is performed on the Basic Calibrated Data 
(BCD) frames, and the FWHM of the point-spread func¬ 
tion is ~1.7" in both bands. 

In addition, the individual frames are mosaicked, and 
p oint sources identified. During catalog construction 
(§ 2.31, the point source lists are also position matched to 
2MASS JHKs photometry. The 2MASS bands have been 
dereddened to account for interstellar extinction. The 
reddening co efficients used for th e LMC can be found 
in Table 1 of Riebel et al. (20121. The photometry for 


Eb-v 

= 0.04 mag ( 

Harris & Zaritsky 2004 

Schlegel 

et al. 

1998 

) and the prescription of 

Glass et a 



tudes using zero points of 6.12 mag (280.9 Jy) for [3.6] 
and 5.63 mag (179.7 Jy) for [4.5] ( Reach et al.||2005] ) . 

The pipeline provides 1 cr photometric uncertainty in 
each band, plotted as a function of source magnitude 


in Figure Throughout this paper, we take these 1 cr 
errors as the uncertainties in our photometry. The dis¬ 
tributions in Fig. show a hard cutoff at 0.2 mag due to 
the requirement that sources in the original SAGE survey 
have a S/N > 5, which effectively results in the top cutoff. 
Sources brighter than magnitude ~9.0 are observed only 
using the 0.6 s exposure mode, which results in a larger 
uncertainty than the 12 s exposure time. Between magni¬ 
tudes 9 and 12, the 0.6 s and 12 s exposures are averaged, 
and the photometry of sources dimmer than magnitude 
^12 is derived only from the 12 s exposures. These three 
regimes of source brightness result in the breaks visible 
at ~9 and ~12 mag. See § 3.1.3 of Meixner et al. (20061 
for a full discussion. 


2.3. Catalog Construction 

The initial data product was a ‘full’ point source list, 
without cosmic ray screening. These artifacts were re¬ 
moved from the catalog by matching the full lists to the 
SAGE mosaic photometry archive^ using a 2" match¬ 
ing radius. The mosaic archives of the original SAGE 
surveys were constructed by co-adding and re-reducing 
the photometry from the original two SAGE Epochs. 
The archives are deeper, more complete, and of higher 
signal-to-noise (S/N) than the SAGE-Var data due to the 
greater exposure time of the mosaic photometry (up to 
^50 s per pixel compared to 12 s) and the fact that the 
original surveys were performed during the cold phase 
of the Spitzer mission whereas the SAGE-Var observa¬ 
tions were taken after the liquid helium cryogen aboard 
Spitzer was exhausted. The higher operating tempera¬ 
ture of the Spitzer instruments restricts our observations 
to only the shortest IRAC wavelengths, and lends addi¬ 
tional thermal noise to our images. Due to their greater 
depth and quality, the SAGE mosaic archives serve as a 
‘truth field’ for the SAGE-Var full source lists, providing 
a more thorough list of actual astronomical sources while 
screening out instrumental artifacts by position match¬ 
ing. After the SAGE-Var source list was matched to the 
mosaic photometry, the original SAGE Epoch 1 and 2 
archive data were matched individually as well, produc¬ 
ing a final source list of 1717 554 (LMC) and 457 760 
(SMC) unique sources with up to 6 individual 3.6 and 
4.5/rm observations, as well as mosaic photometry avail¬ 
able in the longer wavelength IRAC bands. Tablelists 
the number of unique sources detected in each band of 
each epoch of SAGE-Var observations. Also listed is the 
number of objects detected in only one epoch, in two 
epochs, in three epochs, etc. A greater number of de¬ 
tections is preferable. The entire SAGE-Var data set, 
including all 6 epochs of observations in both bands, is 
available as an online table, hosted at the NASA/IPAC 
Infrared Science Archive^. 


2.4. Source Classification 


Variables in the SAGE-Va r dataset were iden tified us¬ 
ing the variability criteria of Vijh et al. (2009). We cal- 


^ http://irsa.ipac.caltech.edu/data/SPITZER/SAGE/doc/ 
\SAGEDataProductsDescriptioii_Sep09 .pdf 

^ http://irsa.ipac.caltech.edu/Missions/spitzer.html 
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Fig. 3.— Each plot shows the la uncertainty produced from the Wisconsin pipeline as a function of source magnitude for the entire 
SAGE-Var dataset as a Hess diagram. The top row shows the data from the LMC ([3.6] on the left, [4.5] on the right), while the bottom 
row shows the SMC data. 


TABLE 2 
Object Count 


Epoch 

LMG 

SMG 


3.6 fim 

4.5 fim 

3.6 fim 

4.5 fim 

Mosaic Photometry 

1,712,135 

1,712,994 

454,951 

455,587 

SAGE Epoch 1 

1,185,774 

1,083,380 

297,732 

269,178 

SAGE Epoch 2 

1,163,888 

1,068,438 

290,055 

262,408 

SAGE-Var 1 

1,177,715 

824,928 

289,359 

185,392 

SAGE-Var 2 

1,242,129 

852,734 

264,921 

176,911 

SAGE-Var 3 

1,244,543 

855,533 

287,168 

173,756 

SAGE-Var 4 

1,244,791 

845,370 

280,656 

172,379 

1 Detection 

116,285 

240,980 

56,602 

78,408 

2 Detections 

197,701 

240,983 

75,528 

61,853 

3 Detections 

270,324 

273,872 

79,241 

65,429 

4 Detections 

253,951 

183,972 

54,222 

45,670 

5 Detections 

271,603 

180,249 

64,538 

43,291 

6 Detections 

593,727 

391,448 

120,822 

73,748 

Total Sources 

1,717 

,554‘‘ 

457,' 

760^ 


Note. — Summary of the number of unique objects detected 
in the SAGE-Var survey of the SMC and LMC 


These numbers represent the total number of unique sources de¬ 
tected in each galaxy, with an absolute minimum of one detection 
in one band. 


culate variability indices 



for every star, for every possible combination of epochs i 
and for the flux / in each SAGE-Var band b, ([3.6] and 
[4.5]) with photometric uncertainty cr. The photomet¬ 
ric uncert ainti es are taken directly from the Wisconsin 
pipeline (§ 2.2). The variability index is thus the number 
of st andard errors by w hich two epochs differ in bright¬ 


ness. 


Vijh et al. (2009) had only the first two epochs of 


the original SaGe survey, but 5 bands of photometry. In 
this study, our situation is reversed, in that we have only 
two bands of photometry, but = 15 possible epochal 
differences. Each epochal difference probes variability on 
a different timescale (Table [^, and we term each epochal 
difference an interval. In order for a source to be clas¬ 
sified as a variable in a given interval, we require it to 
exhibit Scr flux variation in the same direction (that is, 
brightening or dimming) in both the [3.6] and [4.5] bands 
(If^jl — ^ same sign in both bands). 

The SAGE-Var sample h as a total of 819 sources in 
common with the sample of Vijh et al. (2009|). Of these, 
we independently identify 752 (92%) as variables using 
our own criteria. The remaining 67 sources typically 
show marginally variable behavior, with variability in¬ 
dices very close to, but not quite exceeding, our 3(T level. 
We manually add to our catalog 66 of these 67 sources, 
omitting one source only detected in one of the four new 
epochs of SAGE-Var. 

These criteria resulted in 2198 unique variables in the 
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TABLE 3 

Full Object Catalog Contents 


Column 

Name 

Description 

Null 

1 

desig 

Source IRAC Designation 


2 

ra 

Right Ascension, J2000 [deg] 


3 

dec 

Declination, J2000 [deg] 


4 

el_36 

Epoch 1 flux in [3.6] band [Jy] 

-99 

5 

el_36_u 

Uncertainty in Epoch 1 [3.6] flux [Jy] 

-99 

6-15 

eN_36 

Epoch 2-6 flux/uncertainty in [3.6] band^ [Jy] 

-99 

16-27 

eN_45 

Epoch 1-6 flux/uncertainty in [4.5] band'^ [Jy] 

-99 

28 

mean_36 

Mean [3.6] flux [Jy] 

-99 

29 

mean_36-U 

RMS Uncertainty in mean [3.6] flux [Jy] 

-99 

30 

mean_45 

Mean [4.5] flux [Jy] 

-99 

31 

mean_45_u 

RMS Uncertainty in mean [4.5] flux [Jy] 

-99 

32-46 

var_36-N 

Variability index at [3.6] for interval N 

NaN 

47-61 

var_45_N 

Variability index at [4.5] for interval N 

NaN 

62 

ogle_id 

ID From the OGLE-III Catalog of Variable Stars 


63 

ogle.class 

Classification from the OGLE-III CVS 


64 

ogle.per 

Variability period from the OGLE-III CVS [days] 

-99 

65 

macho _id 

ID from the MACHO survey 


66 

macho.per 

Variability period from the MACHO survey [days] 

-99 

67 

grams _class 

Classification of best fitting GRAMS model (G or O) 


68 

y so .class 

‘Y’ if a source is classified as a YSO candidate 



Note. — The full SAGE-Var catalog of 1717 554 (457 760) sources in the LMC (SMC) is 
available from IRSA. This table is provided as a guide to the online catalog’s structure and 
content. 

^ The epoch 2 through 6 [3.6] photometry follows the same format as columns 4 and 5 
^ The [4.5] ph otom etry follows the same format as the [3.6] photometry in columns 4-15 
^ Defined in § |2.4| and Table 


TABLE 4 

Variable Object Catalog Contents 


Golumn 

Name 

Description 

Null 

1 

desig 

Source IRAC Designation 


2 

ra 

Right Ascension, J2000 [deg] 


3 

dec 

Declination, J2000 [deg] 


4 

el_36 

Epoch 1 flux in [3.6] band [Jy] 

-99 

5 

el_36-U 

Uncertainty in Epoch 1 [3.6] flux [Jy] 

-99 

6-15 

eN_36 

Epoch 2-6 flux/uncertainty in [3.6] band^ [Jy] 

-99 

16-27 

eN_45 

Epoch 1-6 flux/uncertainty in [4.5] band'^ [Jy] 

-99 

28 

mean_36 

Mean [3.6] flux [Jy] 

-99 

29 

mean_36_u 

RMS Uncertainty in mean [3.6] flux [Jy] 

-99 

30 

mean_45 

Mean [4.5] flux [Jy] 

-99 

31 

mean_45-U 

RMS Uncertainty in mean [4.5] flux [Jy] 

-99 

32-46 

var_36_N 

Variability index at [3.6] for interval N 

NaN 

47-61 

var_45-N 

Variability index at [4.5] for interval N 

NaN 

62 

ogle_id 

ID From the OGLE-III Catalog of Variable Stars 


63 

ogle.class 

Classification from the OGLE-III CVS 


64 

ogle.per 

Variability period from the OGLE-III CVS [days] 

-99 

65 

macho _id 

ID from the MACHO survey 


66 

macho.per 

Variability period from the MACHO survey [days] 

-99 

67 

grams.class 

Classification of best fitting GRAMS model (C or O) 


68 

y so .class 

‘Y’ if a source is classified as a YSO candidate 


69 

amp_36 

SAGE-Var observed [3.6] Amplitude 


70 

amp_45 

SAGE-Var observed [4.5] Amplitude 



Note. — This table extracts just those sources fiagged as variable by the criteria of § |2.4| 
It follows essentially the same format as Table with the addition of the observed amplitudes 
of the variable objects. These are simply the difference between the brightest and dimmest 
magnitudes observed for the source, and represent a lower limit on the source’s full variability. 
^ The epoch 2 through 6 [3.6] photometry follows the same format as columns 4 and 5 
^ The [4.5] photometry follows the same format as the [3.6] photometry in columns 4-15 
^ Defined in § |2.4| and Table 










6 


Riebel et al. 


TABLE 5 

Interval Timescales 


Epochs 

LMC Interval 
(days) 

SMC Interval 

(days) label 

U 

3 

O 

in 

Epoch 1 — 

Epoch 2 

100 

96 

1 

0) 

Epoch 1 — 

Epoch 3 

1854 

793 

2 


Epoch 1 — 

Epoch 4 

1878 

819 

3 

o 

> 

Epoch 1 — 

Epoch 5 

1984 

922 

4 


Epoch 1 — 

Epoch 6 

2107 

1096 

5 


Epoch 2 — 

Epoch 3 

1754 

697 

6 

0) 

n 

Epoch 2 — 

Epoch 4 

1778 

723 

7 

t 

Epoch 2 — 

Epoch 5 

1884 

826 

8 

z 

Epoch 2 — 

Epoch 6 

2007 

1000 

9 


Epoch 3 — 

Epoch 4 

24 

26 

10 


Epoch 3 — 

Epoch 5 

130 

129 

11 


Epoch 3 — 

Epoch 6 

253 

303 

12 


Epoch 4 — 

Epoch 5 

106 

103 

13 


Epoch 4 — 

Epoch 6 

229 

277 

14 


Epoch 5 — 

Epoch 6 

123 

174 

15 


Note. — 

By taking the difference between all possible com- 

w 

binations of observation epochs, we probe source variability 

u 

on 15 timescales, from 

month to ~ 

5.5 years. Throughout 

3 

o 

this paper, ’ 

we refer to these epochal differences as intervals. 

in 

The “label” 

column refers to the electronic table of the entire 

n 


SAGE-Var catalog, available from the NASA/IPAC Infrared 
Science Archive (IRSA) (columns 32-61, Tables fc [Z]! . 


LMC along with 571 in the SMC. Histograms showing 
the number of sources seen as variable in each interval are 
shown in Figure]^ Intervals 1-9 show a generally higher 
number of detected variables because they all compare 
at least one of the initial two SAGE epochs to later ob¬ 
servations. Taken during the cold Spitzer mission, which 
had a S/N approximately twice that obtained during the 
warm mission when all other epochs in SAGE-Var were 
taken. Interval 10 shows an unusually small number of 
variables because it is an order of magnitude shorter than 
any other interval, spanning only ~20 days. A randomly 
phased interval that short is unlikely to catch variability 
for many infrared sources. 

In order to classify the variables we did detect, we 
matched our detected variables against the OGLE-III 
Catalog of Variable Stars (CVS)^, which consists of 
^150,000 classified variables in the LMC (~26,500 in 
the SMC) with well characterized variability information. 
Again using a 2" matching radius, we find 1361 OGLE- 
CVS matches to our 2198 SAGE-Var variables in the 
LMC, and 323 matches to our 571 variables in the SMC. 
These matches are shown in the [3.6] — [4.5] vs. [4.5] 
CMD in Figure]^ The OGLE populations we find in our 
data are detailed in Tables along with the appro¬ 

priate reference to the relevant OGLE-CVS document, if 
any. See the cited references for thorough discussion and 
definitions of the various classifications. The OGLE ob¬ 
servations were taken primarily in the I band, with some 
additional images in V and B. The OGLE camera has a 
pixel scale of 0.417"/pi xel, and concentrates 80% of the 
light in a 0.5" disk. See|Udalski et al. (1997) for a thor¬ 
ough discussion of the OGLE hardware and observing 
protocols. 

LPVs, mainly evolved AGB stars, are the most numer¬ 
ous OGLE sources we detect in our sample. This is not 
surprising, as AGB stars are among the brightest objects 
in the IR sky, and many of our intervals probe timescales 


http://ogledb.astrouw.edu.pl/-ogle/CVS/ 
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Interval 
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1000 
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Fig. 4.— Histogram of the number of sources detected as vari¬ 
able in each interval of SAGE-Var. The top panel shows the source 
distribution in the LMG, while the bottom panel is that for the 
SMC. The colored portion of each bar represents the number of 
sources classified by the OGLE project. Red: LPVs, Cyan: Clas¬ 
sical Cepheids, Hot Pink: Type H Cepheids, Blue: Eclipsing Bina¬ 
ries, Brown: YSOs, Green: R CrB stars. Purple: RR Lyrae stars. 
The white area at the top of each interval represents the unclassi¬ 
fied variables seen in SAGE-Var but not listed in the OGLE-CVS 
(Figure 1^. The labels for the intervals on the x-axis refer to the 
indices in Table Classes are plotted from least numerous to most 
numerous, bottom to top. 


on which LPVs are expected to vary. Sources identi¬ 
fied as LPVs previously have be en matched to varia bility 
data from the MACHO survey (Riebel et al.||2010|). The 
MACHO survey utilized two non-standard bands; a red 
(^690/tm) and a blue (~520/im). The MACHO optics 
have a pixel scale of 0.635"/pixel and a median PSF of 



fication besides LPVs as well, down to nearly the limiting 
magnitude of our survey. These other classes tend to be 
blue in the diagrams of Figure]^ as the instability strip 
they occupy on the HR diagram places the peak of their 
SEDs blueward of the Spitzer bands. 

We also matched our data to the list of Young Stellar 
Ob ject candidates (YSO s) in the LMC (SMC ) compiled 
by ^ ^ 


50(1 


Carlson et al. 

(H37J 


(|2012|) (iSewilo et al.||2013l. 
YSOs in the entire SAGE-Var LM' 


We hnd 
C (SMC) 


dataset, but only 12 (4) of them are identihed as variables 
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Fig. 5.— [3.6] — [4.5] vs. [4.5] CMD highlighting the OGLE-CVS 
variables with counterparts amongst the SAGE-Var classified vari¬ 
ables. The LMC sources are shown in the top panel and the SMC 
sources are shown in the bottom. The entire SAGE-Var catalog is 
shown as a gray-scale Hess diagram in the background. The same 
color scheme as in Figure]^ is used. 

using our criteria. These are also shown in Figure!^ 
Figure shows the same CMDs as Figure but high¬ 
lighting the sources without counterparts in the OGLE 
CVS or MACHO survey. We remove any source with a 
SAGE [4.5] mosaic photometry dimmer than 15**' mag¬ 
nitude identified as a variable but without an OGLE or 
MACHO counterpart. We include these sources in the 
online catalog for completeness, but data artifacts (such 
as blending with nearby sources) cause problems with our 
variability criteria for such dim sources. After removing 
those dim sources, we are left with 641 (139) IR variables 
in the LMC (SMC) without OGLE or MACHO detected 
variation. There are a few extremely red sources in the 
AGB region of the CMD without OGLE identifications. 
Most of these sources are flagged as probable variables by 
the WISE survey, and our results bolster that conclusion 
(§| 0 . 


3. RESULTS 
3.1. New AGBs 


Fig. 6.— [3.6] — [4.5] vs. [4.5] GMD highlighting variables de¬ 
tected by SAGE-Var but not identified by the OGLE or MACHO 
surveys. The LMC sources are shown in the top panel, and the 
SMC sources are shown in the bottom. The entire SAGE-Var cat¬ 
alog is shown as a gray-scale Hess diagram in the background. To 
eliminate false positives due to low S/N, we only consider sources 
brighter than magnitude to be variables. 


treme AGBs, exactly the stars whose variability is most 


dominate the mass return from AGBs to the ISM ( 

Riebel 

et al. 

2012 

Boyer et al. 

2012 Matsuura et al. 

2009). 

SAGE-Var represents the hrst large-scale variability sur- 


vey at such red wavelengths. As such, we have the ability 
to detect variable stars which have never been catego¬ 
rized as such before. Our survey detects variability in 
641 (139) sources in the LMC (SMC) which do not have 
well-defined variability measurements in the OGLE or 
MACHO surveys (Fig. [^. Out of these newly identified 
variables, 10 sources in the LMC and 6 in the SMC can 
be classihed as AGB c andidates, based on previous stud¬ 
ies (Riebel et al. 20121), or their position in the J — Kg vs. 


ATs CMD. Figure!^ 
didates against a 
sample. 


ighlights these 10 (6) new AGB can- 
ess diagram of the entire SAGE-Var 


3.1.1. LMC AGB Candidates 


Because they typically require large amounts of tele¬ 
scope time and the expense of space-based multi-epoch 
observations often makes such lengthy surveys cost- 
prohibitive, most variability surveys are conducted from 
the ground, and hence at optical wavelengths. Ground- 
based variability surveys thus miss the reddest, most ex- 


We find 10 AGB candidates in the LMC without previ¬ 
ously well-measured periods in either the OGLE or MA¬ 
CHO surveys. Most of these were flagged as highly likely 
to be true variables by the WISE survey, and we confirm 
that measurement. We do not find any of these 10 can¬ 
didates in the SAGE-Spec ([Kemper et al.j 2010|) list of 
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TABLE 6 

Matched variable populations detected in SAGE-Var LMC 


Variable Classification Number 


Reference/Definition 


Total Variables Detected 

2198 

Long-Period Variables 

1065 

AGB C 

5 

AGB O 

6 

Mira AGB C 

426 

Mira AGB O 

143 

OSARG AGB G 

15 

OSARG AGB O 

42 

OSARG RGB C 

1 

OSARG RGB O 

17 

SRV AGB G 

262 

SRV AGB O 

148 

Cepheids 

28 

to 

1 

10/20 

1 

F 

26 

Type II Cepheids 

19 

BLHer 

1 

RVTau 

12 

WVir 

6 

RR Lyrae 

3 

R GrB 

6 

Eclipsing Binaries 

25 

EC 

8 

ECL 

1 

ED 

8 

ESD 

8 

Young Stellar Objects 

12 


[Soszynski et al.] l |2009b[ | 


OGLE small-amplitude variable reg giant (OSARG) 


Semi-regular variable (SRV) 

|Soszynski et al.| (|2008a|| 

overtone pulsation mode 

Mixture of and 2^“^ overtone pulsation modes 
Fundamental-mode puls ation 
[Soszynski et al.] l |2008b[ | 


Soszynski et al. 

2009a 

Soszynski et ai. 

1 2uuyc 

UiraczyR et al.| l 

iuLljl 


Ellipsoidal contact binaries 
Detached binaries 
Semi-detached binar ies 
l|2012|l 


Carlson et al. 


Note. — This table lists the variable c lassifications and sub classes used by the OGLE Catalog 
of Variable Stars and the YSO catalog of|Carlson et al.|^2012|, and the numbers of such sources 
also detected as variable by the SAGE-Var survey in the LMC. A “I" matching radius was used for 
all OGLE catalog comparisons, and the YSO list was matched based on SAGE Archive designa¬ 
tions. See the cited references for complete definitions of the subcategories listed here. Uncommon 
abbreviations are defined in the third column. 


TABLE 7 

Matched variable populations detected in SAGE-Var 
SMC 


Variable Classification 

Number 

Reference 

Total Variables Detected 

571 


Long-Period Variables 

276 

[Soszynski et al.j ([2011b[) 

Mira C 

140 


Mira O 

22 


OSARG C 

7 


OSARG O 

4 


SRV C 

95 


SRV O 

8 


Cepheids 

F 

42 

40 

[Soszynski et al.j ([2010a[l 

lO 

2 


Type II Cepheids 

3 

[Soszynski et al.j ([2010b[) 

HV'i'au 

1 


WVir 

1 


pWVir 

1 


Young Stellar Objects 

4 

Sewilo et al. (2013[l 


Note. — This table lists the variable classifications and 


subclassifications used by th e OGLE Catalog of Variable 
Stars aud the YSO catalog of Sewilo et al. (20131, and the 
numbers of such sources also detected as variable by the 
SAGE-Var survey in the SMC. A 2" matching radius was 
used for all OGLE catalog comparisons, and the YSO list 
was matched based on SAGE Archive designation. 


LMC sources. Without a chemical classification derived 
from observed spectra, we therefore employ an Spectral 
Energy Density (SED)-based chemical classification. 

We fit the SEDs of these sources with synthetic SEDs 


from the GRAMS model grid (ISargent et al. 120111 ISrini- 

. 1 V A -I. m-r^ . II I 'l 


vasan et al.|2011 ). GRAMS is a grid of dusty AGR/RSG 


star m odels corhputed by t he radiative transfer code 
2Dust dUeta & Meixner|2003|) . The grid is constructed by 
placing spherically symmetric circumstellar dust shells of 
varying inner radii, optical depth, and both 0-rich and 
C-rich chemical composition around model stellar pho¬ 
tospheres. The models assume a constant mass-loss rate 
(be., an inverse-square density law) in the shell and a 
power-law distribution of grain sizes with an exponen¬ 
tial falloff, with a typical size of ~0.1/im. The 0- rich 
dust is modeled usi ng the astronomical silicates of |Os-| 


senkopf et al. (19921, and the carbonaceous dust is a 9:1 
mixt ure by mass of amo rphous carbon (optical constants 
from Zubko et al.||199 6|) and silicon carbide (optical con¬ 
stants from Pegourie p988 ). We perform minimum 
fitting to obtain the best-fit parameters, which include 
the luminosity, dust-production rate and chemical type. 

Using the GRAMS model grid, we classify 8/10 of these 
AGB candidates as C-rich. C-rich AGB stars tend to 
be redder than 0-rich stars, and our newly measured 
variables skew red. AGB stars in our sample detected 
by the OGLE survey have an average [3.6] — [4.5] color 
of —0.07 mag, while the sources with no OGLE detection 
are nearly a magnitude redder at these wavelengths, with 
an average [3.6] — [4.5] of 0.74 mag. SAGE-Var has too 
sparse a sampling of the lightcurves to deduce a period 
for their variability, but we can place lower limits on 
their IR variability amplitudes. These sources are listed 
in Table m 
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TABLE 8 

AGB Candidates in the LMC without OGLE or MACHO Variability Measurement 


SAGE 

Designation®' 

RA 

(2000) 

Dec 

(2000) 

GRAMS 

Class 

[3.6] 

Amplitude 

[4.5] 

Amplitude 

[3.6]-[4.5] 

color 

J051041.21-683606.6 

77.6717 

-68.6018 

c 

0.83 

0.85 

1.24 

J051414.85-700409.8 

78.5619 

-70.0694 

c 

0.93 

0.86 

1.12 

J052503.26-692617.3 

81.2636 

-69.4381 

c 

0.90 

0.69 

1.48 

J052813.02-691228.4 

82.0543 

-69.2079 

c 

0.86 

0.75 

1.09 

J052900.19-695247.3 

82.2508 

-69.8798 

c 

1.11 

0.52 

0.80 

J053051.75-694328.0 

82.7156 

-69.7245 

c 

1.42 

1.19 

1.20 

J050202.38-690726.2 

75.5099 

-69.1239 

c 

0.52 

0.56 

1.10 

J050718.89-683850.4 

76.8287 

-68.6474 

c 

1.50 

1.27 

1.61 

J051913.89-693818.3 

79.8079 

-69.6384 

o 

0.29 

0.31 

0.26 

J053010.30-690933.8 

82.5429 

-69.1594 

o 

0.36 

0.30 

1.05 


Note. — The 10 AGB candidates in the LMC with variability newly detected by 
SAGE-Var. Many of these sources are confirmed to be AGB stars in the literature, 
but none have previously been observed to vary. The GRAMS Class column lists 
the classification (0-rich or C-rich) assigned each source by the GRAMS model grid. 
The next two columns represent a lower bound on the IR variability amplitude of the 
sources, the maximum observed magnitude minus the minimum observed magnitude. 
The [3.6] — [4.5] color is also listed to connect the entries in this table to Figure]^ 
Designations in the online data table are prefaced with ‘SSTISAGEMA’ 
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Fig. 7.— [3.6] — [4.5] vs. [4.5] CMD highlighting newly identified 
LPV candidates. The LMC sources are shown in the top panel, and 
the SMC sources are shown in the bottom. The entire SAGE-Var 
catalog is shown as a gray-scale Hess diagram in the background. 
Stars classified as C-rich by the GRAMS model grid are shown in 
red, stars classified as O-rich are shown in blue. 


3.1.2. SMC A GB Candidates 

After cross-matching our list of AGB candidates in the 
SMC with SIMBAD and removing all sources confirmed 
to not be AGB stars, we generate a list of 6 new AGB 
candidates. We verified that none of these candidates 
were in the Ruffle et al. (2015, in prep.) list of spec¬ 
troscopically classified SMC sources. Using the GRAMS 
model grid, we classify 5/6 of them as C-rich. These 
sources are listed in Table [9l Source SSTISAGEMA 
J010041.61-723800.7 is the reddest new LPV candidate 
we idenitfy in the SMC. The best-fitting GRAMS model 
is an O-rich model, but the extreme redness of the source 
is more consistent with a C-rich star. Currently, the 
GRAMS model grid in the SMC places too much weight 
on the 24 /rm photometry, and this leads to some sources 
being mis-classified (see Srinivasan, et al. 2015, in prep 
for details). We manually change the classification of this 
source to C-rich. 


3.2. Variability Amplitude and Dust Production Rate 


Le Bertre (1992) found a strong correlation between 
variability ainplitude (measured in the K band), and 
K — L' col or for a sample of 20 carb on-rich LPVs in 
our galaxy. iWhitelock et al.|(|1991|) and|Le Be rtre] (|1993|) 

rtor U-rich LPVs. IWhitefocH 


found a siniilar correlation 


et al. (1994) found a strong correlation between the IR 
amplitude and the K — [12] color. The K band is dom¬ 
inated by light from the stellar photosphere, while the 
IRAS 12 /rm band represents the emission from the cooler 
dust shell around the star, and so this color serves as an 
indicator of the thickness of the circumstellar dust shell, 
and hence the rate of dust production by the star, if dust 

production rate (DPR) is assumed to be const ant. _ 

By using the GRAMS modeling results of|Riebel et 


a y usin g tne GKAivib moaeimg results oi iKiepei et 
al. i|2012 ), we can directly compare the [3.6] amplitude 


and the DPR for our sample. Our results are shown in 
Figure [8] We extract every source in the SAGE-Var cat¬ 
alog with a valid GRAMS cl assifi cation which was also 
classified as a true variable (§) 
relation of increasing DPR wit 


>2.4). We see a slight cor- 
;h increasing infrared am¬ 
plitude (albeit one with a very large scatter). G-rich 
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TABLE 9 

AGB Candidates in the SMC without OCLE or MACHO Variability Measurement 


SAGE 

Designation®' 

RA 

(2000) 

Dec 

(2000) 

GRAMS 

Class 

[3.6] 

Amplitude 

[4.5] 

Amplitude 

[3.6]-[4.5] 

color 

J005106.28-731635.9 

12.7762 

-73.2767 

G 

0.34 

0.34 

0.14 

J004544.12-720815.4 

11.4338 

-72.1376 

G 

0.32 

0.40 

0.24 

J005926.35-722341.4 

14.8598 

-72.3949 

G 

1.09 

0.88 

1.23 

J010232.75-721912.5 

15.6365 

-72.3202 

G 

1.08 

1.06 

1.17 

J010041.61-723800.7 

15.1734 

-72.6335 

O'" 

0.95 

1.03 

1.69 

J005131.21-732007.7 

12.8801 

-73.3355 

O 

0.34 

0.23 

-0.06 


Note. — The 6 AGB candidates in the SMC with variability newly detected by 
SAGE-Var. Most of these sources are confirmed AGB stars, but none have previously 
been observed to vary. The GRAMS Class column lists the classification (0-rich or C- 
rich) assigned each source by the GRAMS model grid (Srinivasan, et ah, 2015, in prep). 
The next two columns represent a lower bound on the IR variability amplitude of the 
sources, the maximum observed magnitude minus the minimum observed amplitude. 
The [3.6] — [4.5] color is also listed to connect the entries in this table to Figure]^ 
Designations in the online data table are prefaced with ‘SSTISAGEMA’ 

The best-fitting GRAMS model for this source is an O-rich model with poor fit 
quality. Based on the extreme redness of this source, we classify it as C-rich 


10 "' 

10 "° 

10 "° 

10 "'° 

10 "" 

10 "'° 

0.0 0.5 1.0 1.5 2.0 

[3.6] Amplitude 

Fig. 8.— Dust Production Rate (DPR, Mq yr“^) vs [3.6] Am¬ 
plitude for AGB stars in the LMC and SMC. The DPR is taken 
from the best-fitting GRAMS model, and the [3.6] Amplitude is 
defined as the dimmest measured SAGE-Var [3.6] magnitude mi¬ 
nus the brightest [3.6] magnitude. C-rich AGB stars are shown as 
open red circles, and O-rich AGB stars as filled black triangles. 


DPR VS [3.6J Amplitude 



AO-rich AGB stars 
O C-rich AGB stars 


and O-rich AGB stars have considerable overlap in their 
range of DPR, but more C-rich stars tha n O-rich stars 
exten d to higher DPRs (see discussion in Riebel et al. 
20121. Because of this, the trend is more visible among 
the C-rich sources, but both populations follow the same 
basic trend, with considerably more scatter at low am¬ 
plitudes for the O-rich sources. 

3.3. Period-Luminosity Relationships 
3.3.1. AGB Period-Luminosity Relationships 

Characterizing the AGB Period-Luminosity (PL) rela¬ 
tionship is important in light of the upcoming launch of 
the James Webb Space Telescope (JWST). AGB stars 
are among the brightest sources in the infrared (IR) sky, 
and with a large orbital platform concentrating on the 
IR, AGB stars could serve as important distance indi¬ 
cators, if their intrinsic luminosity can be well deter¬ 
mined from their variability. Previous investig ations of 
the AGB PL relationship (see Riebel et al.|2010 and ref¬ 
erences therein) have found a scatter of a few tenths of 
a magnitude about the best-fit line, and it has not been 


clear whether this scatter was intrinsic to the relationship 
(i.e. a real phenomenon due to an astrophysical cause 
such as a spread in stellar masses) or an artifact result¬ 
ing from observing stars at various points in their light 
curves. Constructing mean magnitudes from observa¬ 
tions taken at several points in an AGB star’s lightcurve 
are an essential step toward clarifying the nature of this 
spread. The 6 epochs of the SAGE-Var survey allow us 
to calculate improved mean magnitudes for every source 
in the catalog. 


We use only AGB stars from the sample of Riebel et 


al. (2010) which are identified as variables in the SAGE- 


Var data to construct our relations. Only stars that have 
been observed at more than one brightness benefit from 
the averaging process, as opposed to stars with 6 mea¬ 
surements of the same brightness. For all AGB stars 
identified as variables, all available valid flux measure¬ 
ments from the SAGE and SAGE-Var surveys were av¬ 
eraged to construct mean fluxes at both [3.6] and [4.5]. 
Photometric errors from each observation were added in 
quadrature to produce an error for the mean flux. These 
mean fluxes were tra nsfor med into magnitudes using the 
zero points given in § |2.2| 

Figure [^ shows the PL relations constructed from the 
6-epoch mean magnitudes constructed from the SAGE- 
Var data. The 3.6 fim magnitude is used in the top panel, 
and the 4.5 fim band is used in the bottom. The periods 
come from the MACHO survey if possible, and from the 
OGLE survey if no MACHO period was available. The 
MACHO periods were given pr iority in order to ma intain 


I 


For 

GLE 


consistency with the study of Riebel et al. (2010 
those stars with measured periods m both the 
and MACHO surveys, the periods were found to differ 
by only 3 days on average. We only present updated fits 
for stars pulsating in the fundamental mode (sequence 
1), as this was the only sequence with a significant pop¬ 
ulation observed in SAGE-Var. The derived PL rela¬ 
tions for the stars classified as 0-riclMby the best-fitting 
GRAMS model) are shown in FigurelWand quantitatively 
described in Table [TOl The fits to the stars identified as 
C-rich are shown in Figure [TO] and described in Table [TT| 
In Figurej^stars with periods > 300 days are systemat- 





















TABLE 10 

Period Magnitude Relationships eor 0-rich AGB stars 
Pulsating in the Fundamental Mode 
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PL Relation 

N 

scatter (mag) 

This work: 

[3.6] = -4.00(±0.03) log P + 20.24(±0.07) 
[4.51 = -3.78 ±0.03') loe P + 19.72f±0.06') 

131 

131 

0.18 

0.22 

|Riebel 

et aL| 1 2010||: 

2221 

2227 

0.271 

0.270 

[j-o] = 

[4.5] = 

-d.4l(±u.u4) log P ± 18.88(±0.09) 
-3.35(±0.04) log P + 18.80{±0.01) 


Note. — PL relations for LPVs classified as O-rich AGB 
stars. The quoted scatter is the standard deviation of the 
residuals about the best-fit line. 


ically brighter than the trend line, whereas in Figure [10| 
stars with periods < 300 days are systematically dimmer 
than the trend line. These outlier stars are an overlap 
between the two populations. The C-rich best-fit trend 
line is about a magnitude brighter than the O-rich fit, 
and C-rich stars have longer periods on average. Thus, 
the shortest period C-rich stars that are dimmer than the 
trend line actually overlap the main body of the O-rich 
stars in PL space. In the same manner, the longest pe¬ 
riod O-rich stars which are brighter than the trend line 
in figure o verlap the main population of C-rich stars 
in figure [T^ When using photometric classification, the 
division between O- and C-rich AGB stars is not sharp, 
and the two populations smoothl y blend into one an other 
in PL space (see Figures 1 & 5 in Riebel et al.|20r0 ). Be¬ 
cause there are only about 10 of these outlying stars, the 
fits are insensitive to their inclusion, and they have not 

been removed from the sample prior to fitting. _ 

Using the criteria described in Appendix A of |RiebeI| 


et al. (2010), we identify three O-rich AGB stars in the 
biVlC which lie on Sequence 1. We corrected for the dif¬ 
ference in distance between the LMC and the SM C using 


distance modu li of 18.54 and 18.93, respectively (Keller 


& Wood 2006). These stars, shifted to the distance of 


the LMC, are plotted as pink crosses in Figure On 
average, the LMC has a greater metallicity than that of 
the SMC. Our numbers are too small to draw definite 
conclusions, but we do not see any evidence for a de¬ 
pendence of the AGB PL relation on metallicity in our 
sample. This idea is worthy of further investigation. 

While we do see a slight decrease in the scatter about 
the be st-fit line compared to that found by [Riebel et al.| 
(2010), the reduction is less than a factor of two, which 
leads us to believe much of this remaining scatter is in¬ 
trinsic to the relationship and is not a product of observ¬ 
ing different stars at different phases of their lightcurve. 

We identified 223 (49) Carbon-rich, non-extreme AGB 
stars pulsating in the fundamental mode in the LMC 
(SMC). “Extreme” AGB stars are defined as those stars 
with J — [3.6] > 3.1. The flux in the J band traces 
the emission from the stellar photosphere, while the [3.6] 
emission is from the circumstellar dust shell around these 
heavily enshrouded stars. At SAGE-Var wavelengths, 
these stars do not follow the same PL relation as less en¬ 
shrouded stars because their brightness does not reflect 
stellar emi ssion, but emissio n from the dust shell (see 
figure 4 in Riebel et al.ll2010| . Sources thus classified as 
extreme are not incfuded in this fit. The derived PL re la- 
tions to the remaining G-rich stars are shown in Figure [Tol 



PL Relation for 0 —rich LPVs at [4.5] 



Period (days) 

Fig. 9. — PL relationships for GRAMS classified O-rich AGB 
stars in the SAGE-Var sample, constructed using 6-epoch mean 
magnitudes in the [3.6] and [4.5] bands. Stars from the LMG are 
shown as black points, stars from the SMG shown as pink plus 
signs. Note the horizontal axis is a log scale. 


and quantitatively described in Table 11 The [4.5] data 
are shown for completeness, but we do not attribute any 
sig nificance to the line ar fits in that band. As discussed 
by [Riebel et al.| (|2010[), AGB stars are affected by a CO 
absorption feature m the [4.5] band which distorts the 
linear PL relationship in that band. O-rich stars do not 
have enough CO in their atmospheres t o be greatly af¬ 
fected. Figure 6 in Riebel et al. (2010) shows that the 
slope of the O-rich AgB PL relation is relatively insen¬ 
sitive to wavelength, whereas the C-rich relation shows 
significant variation. 

Correcting for distance, the SMC stars are well-mixed 
with the LMC stars, showing comparable scatter about 
the best-fit line. With many more samples than the O- 
rich, this provides more substantial evidence that C-rich 
AGB stars in the SMC follow the same PL relationship 
as those in the LMC. 


3.3.2. Cepheid IR PL Relationships 

The well-known period-luminosity relationship of clas¬ 
sical Cepheid variables (the Leavitt Law) is one of the 
most important extragalactic distance indicators in as¬ 
tronomy. C epheid studies have gen erally focused on the 
optical, but McGonegal et al. (1982) pointed out that the 
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TABLE 11 

Period Luminosity Relationships for C-rich AGB stars 
Pulsating in the Fundamental Mode 


PL Relation 


N scatter (mag) 


This work: 

[3.6] = -3.63(±0.02)logP + 19.11{±0.05) 272 

[4.5] = -3.84(ifc0.0 2] log P + 19.66(ifc0.06] 272 

iRiebel et al.l ( 20Tol: 

[3.6] = -3.V7(±U.i/5) log P + 19.35(±0.12) 1813 

[4.5] = -3.56(±0.05) log P + 18.96{±0.12) 1816 


0.24 

0.30 

0.251 

0.265 


Note. — PL relations for LPVs classified as C-rich AGB 
stars. The quoted scatter is the standard deviation of the 
residuals about the best-fit line. 


PL Relation for C-rich LPVs at [3.6] 



Period (days) 


PL Relation for C —rich LPVs at [4.5] 



Period (days) 

Fig. 10.— PL relationships for GRAMS classified G-rich AGB 
stars in the SAGE-Var sample, constructed using 6-epoch mean 
magnitudes in the [3.6] band. Stars from the LMG are shown as 
black points, stars from the SMG shown as pink plus signs. Note 
the horizontal axis is a log scale. 


near-IR offers several advantages for calibration of the 
Cepheid PL relation. Interstellar reddening, which im¬ 
poses an intrinsic scatter on magnitude determinations in 
the optical, is much smaller at Spitzer wavelengths. Also, 
Cepheid variability amplitudes are considerably smaller 
in the IR, and thus sparsely sampled lighcurves are good 
tracers of the mean brightness. For these reasons, several 


PL relation in the IR i 
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Ngeow et a 
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Cepheids, with well sampled light curves. SAGE-Var was 


designed as a more general variability survey, so what we 
lack in thorough coverage of the light-curve, we make up 
for in increased sample size. Based on OGLE classifica¬ 
tion, we find 837 (1536) fundamental mode Cepheids in 
our LMC (SMC) data. 

We extract every source observed by SAGE-Var and 
classified as a fundamental mode Cepheid by the OGLE 
survey (Cepheids primarily varying in the first or higher 
overtone mode will by definition not lie on the primary 
PL relation). We eliminate any source with only one 
valid observation during SAGE-Var, as this provides lim¬ 
ited means to constrain the uncertainty on the mean flux 
of the source. We calculate a simple mean magnitude for 
all our sources by averaging the fluxes of all SAGE-Var 
observations and then converting to magnitudes. Only 4 
of our sources are in common with Scowcroft et al. (2011), 
and our mean magnitudes derived from randomlyfphased 
observations are within 0.05 magnitude s of theirs. We 


use the period s from the OGLE survey (Soszyhski et al. 


2008a 2010a). 

'The least squares linear fits (using 3cr clipping) are 
given below in Table [1^ The Leavitt Law at [ 3.6] derived 
from the LMC sampleis illustrated in Figure [TT] and the 
relation at [4.5] is shown in Figure 12 Fits are of the 
form y = A log P + B, with y the mean magnitude of the 
source, and the period P measured in days. The scat¬ 
ter is defined as the standard deviation of the residuals 
to the fit. We visually inspected all the sources clipped 
as part of the fitting process. More than 90% are obvi¬ 
ous blending/confusion issues in the original data. We 
find that when fitting the SMC data, using only sources 
with logP >0.5 provided a bett er visual fit to the dat a 
(this method was also used by Scowcroft et al. 2011). 
The LMC fits were robust to this decision, and no pe¬ 
riod s elec tion criteria were applied. The results given in 
Table [^ show the greater scatter about the SMC rela¬ 
tionship compared with that in the LMC. This is consis¬ 
tent with the much greater overall line-of-sight depth of 
the SM C compared with LMC seen in the early Cepheid 
work of Caldwell & Coulson (1986), as well as in other 
stellar components. 'This is probably due to gravitational 


al. 

2015 


['he s 


and references therein). 


ope of our [3.6] re lations show very good agree¬ 
ment with t hose found by Freed man et al. (2012) in the 
LMC and by Ngeow et al. (2012) in the SMC. I n the [4.5] 
band, we ov erplot the relations determined by|Scowcroft| 


et al. (2011), which also show good agreement with our 

I ' 

bi 


own. As discus sed by [Scowcroft et al. ([2011) and Mon- 
son et al. (20121, CO absorption in the [4.5] band renders 


this PL relation problematic for Cepheid distance deter¬ 
minations. This effect is clearly seen in Table where 
the slopes of the PL relations at [3.6] in the me LMC 
and SMC agree well, while the slopes at [4.5] do not. 
This seems likely to be due to the different effects of CO 
absorption at different metallicities. 


4. CONCLUSIONS 

We present the results from a 4-epoch unbiased IR sur¬ 
vey of the central regions of the LMC and SMC. We have 
produced full catalogs of our observations, consisting of 
1717 554 (457 760) objects in the LMC (SMC). We have 
identified 2198 (571) objects in the LMC (SMC) as prob¬ 
able variables. 
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TABLE 12 

Classical Cepheid Period Luminosity relationships at [3.6] and 
[4.5] IN the Magellanic Clouds 


PL Relation 


N scatter (mag) 


LMC: 

[3.6] = -3.271(±0.004)logP + 

[4.5] = -3.157(±0.004)logP + 
SMC: 

[3.6] = -3.261(±0.006)logP + 
[4.5] = -3.437(±0.004)logP + 


15.993(±0.003) 

811 

0.13 

15.877(±0.002) 

820 

0.13 

16.511(±0.004) 

452 

0.18 

16.665(±0.002) 

454 

0.18 


Note. — PL relations for fundamental mode classical Cepheids in the 
LMC and SMC. Fits were determined using a standard 3cr clipping pro¬ 
cedure. The fits to the SMC data were determined by only considering 
Cepheids with log P > 0.5. The scatter about the fit is defined to be the 
standard deviation of the residuals. 


[4.5J Cepheid PL Relation in the LMC 



log(P/day) 


[3.6J Cepheid PL Relation in the LMC 



log(P/doy) 


[3.6J Cepheid PL Relotion in the SMC 



log(P/doy) 


Fig. 11.— 3.6/im Leavitt Law for 811 (452) classical Cepheids in 
the LMC (SMC). The LMC is shown in the top panel and the S MC 
in the bottom. The fits are quantitatively detailed in Table |12| 
Stars shown in red have residuals to the fit greater than 3cr ana 
did not contribute to the fit. The relation derived in this study is 
pl otted in blue. For the LMC (SMC) data, th e relation determined 
bylFreedman et al.|(|2012|| lINgeow et al.|2012|l is overplotted in red 
as a dashed ime. 


We identify 10 (6) variable AGB candidates in the 
LMC (SMC) without well- dete rmined variable periods 
from OCLE or MACHO (§3.1). Most of these sources 
were flagged as potential variables in the WISE survey, 
and our independent measurement confirms that proba¬ 
bility. 

Using mean magnitudes constructed from 6 epochs of 
observations, we investigate the PL relationship for ACB 


[4.5J Cepheid PL Relation in the SMC 



log(P/day) 


Fig. 12.— 4.5 /im Leavitt Law for 811 (452) classical Cepheids in 
the LMC (SMC). The LMC is shown in the top panel and the SMC 
in the bottom. The fits are quantitatively detailed in Table |12| 
Stars shown in red have residuals to the fit greater than 3a and 
did not contribute to the fit. The relation deri ved in this study 
i s plot ted in blue. The relation determined by [Scowcroft et al.| 
||2011|l using the original SAGE LMC data is overplotted m red as 
a dashed line. 


stars pulsating in the fundamental mode. We find no 
significant reduction in the scatter about the bes t -fit re - 
lation compared to the results of Riebel et al. (2010), 
indicating the sc atter (^0.2 mag) might be intrinsic to 
the data (§ 
tionships o 


3.3. iL We find no evidence that the PL rela- 
AcJl stars in the SMC (Z~0.04) and in the 
LMC (Z~0.08) are different due to metallicity effects. 

We present infrared PL relations for a sample of 811 
(452) Cepheids in the LMC (SMC). Cepheid amplitudes 
are small in the IR, and our random-cadence result s 
compare very well wit h th ose ofiScowcroft et a l. ( 2011 ), 
Freedman et al. (2012) and IN geowet al. l 2012[ ) (§3.3.2). 
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